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Abstract 

We show that parton confinement in the final state generates large 1/Q^ cor¬ 
rections to Bjorken scaling, thus leaving less room for the logarithmic correc¬ 
tions. In particular, the x-scaling violations at large x are entirely described 
in terms of power corrections. For treatment of these non-perturbative effects, 
we derive a new expansion in powers of for the structure function that 
is free of infra-red singularities and which reduces corrections to the leading 
term. The leading term represents scattering from an off-mass-shell parton, 
which keeps the same virtual mass in the final state. It is found that this 
quasi-free term is a function of a new variable x, which coincides with the 
Bjorken variable x for ^ oo. The two variables are very different, how¬ 
ever, at finite . In particular, the variable x depends on the invariant mass 
of the spectator particles. Analysis of the data at large x shows excellent 
scaling in the variable x, and determines the value of the diquark mass to be 
close to zero, x-scaling allows us to extract the structure function near the 
elastic threshold. It is found to behave as F 2 (1 — Predictions for the 

structure functions based on x-scaling are made. 
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I. INTRODUCTION 


Consider the inclusive scattering of a high energy electron, e + N —>■ e' + X, from a 
nucleon of mass M. The spacelike 4-momentum transferred to the target is g = (z/, q). For 
an unpolarized target the double-differential cross-section is determined by two structure 
functions v) and hF 2 (Q^, (or Fi = MWi, F 2 = ZZIF 2 ), where These 

structure functions are given by the imaginary part of the forward Compton amplitude of 
the virtual photon with 4-momentum q (Fig. |^), where the nucleon vertex F is shown in 
Fig. The hrst diagram in Fig. corresponds to the Impulse Approximation (lA), and the 
second diagram describes the Final State Interaction (FSI) of a struck quark with spectator 
quarks and gluons (these are shown explicitly in Fig. ^). The lA term is expected to become 
dominant in the structure function for > 10 (GeV/c)^. As a result Fi and F 2 turn out 
to be functions of the Bjorken variable x = Q^/2Mz/, i.e. Fi[x, Q"^) —> Fi{x), where Fi{x) is 
directly related to the parton distribution q(a:) (the Bjorken scaling), as for instance 

= F2{x, ^ e‘^xqi(x). ( 1 . 1 ) 

i 

The existing data in fact, show considerable Q^-dependence of the structure functions, as 
for the proton structure function from the BCDMS and SLAC [Q experiments shown in 
Fig. 1^. (The solid lines correspond to a 15 parameter £t [^). Usually the scaling deviation 
of the structure functions is attributed mainly to the logarithmic corrections from gluon 
radiation, which is contained in the first diagram in Fig. Hi- Corrections ~ l/Q^ arising 
from the second diagram in Fig. |I| (higher twist terms) are usually considered as playing a 
minor role in scaling violation, even at moderate Q^. 

This common disregard of the FSI terms for ^ 10 (GeV/c)^ is very surprising, 
especially in view of parton conhnement. At hrst sight, the conhning interaction of partons 
in the hnal state should inhuence the structure function strongly. Consider for instance 
the example of two nonrelativistic “quarks” of mass m interacting via a harmonic oscillator 
potential [Q. These quarks are never free and therefore the system in the hnal state possesses 
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a discrete spectrum. As a result the structure function, F{q, z/), as a function of z/, is given 
by a sum of 5 - functions. Obviously, it looks very different from the structure function 
obtained in the lA, which considers the struck parton as a free particle in the hnal state. 
This paradox can be resolved by introducing a scaling variable y 0 



( 1 . 2 ) 


where (m — ?/)/M is a non-relativistic analogue of the Bjorken variable x. Then expanding 
the structure function F{q,y) = F{q,i>) in powers of 1/g, one hnds in the limit g —> cx) 
and y =const that it becomes a smooth curve, F{q, y) ^ 0 ( 2 /), which coincides with a free 
parton response |^,|^. Although this result appears to conhrm the parton model picture, 
it does not imply that the interaction in the hnal state is not important. The latter has 
been merely incorporated in FQ{y) by an appropriate choice of the scaling variable ?/, which 
diminishes the contribution from higher-order (~ 1/g) correction terms. For instance, a 
different choice of the scaling variable could result in very large or even singular corrections 
to the structure function. 

It is practically impossible to calculate the structure function including FSI, except 
for a few simple nonrelativistic models [^0. Therefore it is very important to hnd an 
optimal expansion of the structure function of conhned systems, which reduces higher order 
corrections (~ 1/g) to the zero order (quasi-free) term. For a non-relativistic case such 
an expansion, which leads to the scaling variable y, Eq. (E3). was proposed by Gersch, 
Rodriguez and Smith p. (This expansion was designed for weakly bound systems, but 
appears in fact to be applicable to conhned systems as well [0). Unfortunately, the situation 
in the relativistic case looks very different, and no simple extrapolation of the non-relativistic 
results seems to be possible. 

In this paper we attempt an optimal expansion for the relativistic structure function, that 
can be applied to conhned systems, i.e., it is free of infra-red singularites, and diminishes 
the contribution from FSI. Then the zero order (quasi-free) term ehectively incorporates 
ehects of confining FSI, and thus can be considered a good approximation for the structure 
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function, valid also for nonasymptotic Q^. Such a quasi-free approximation leads to scaling 
of the structure function in a new scaling variable. Finally, we perform the analysis of data 
in terms of this scaling variable and compare the results with a standard approach. 

The plan is as follows: A relativistic expansion of the structure function in powers of 
1/Q^ is derived in Section 2. It is shown that infra-red singularities generated by conhning 
FSI are eliminated. The leading, quasi-free term is discussed in Section 3. We demonstrate 
there that the structure function is a function of a new scaling variable. The evaluation of 
the hrst correction term for the linear-rising potential is presented in Section 4. An analysis 
of data and predictions for new experiments are given in Section 5. The last Section is 
summary. 


II. RELATIVISTIC STRUCTURE FUNCTION 

Consider the nucleon structure function W given by the imaginary part of the forward 
Compton amplitude. Fig. |^. For the sake of simplicity we take all the consituents and the 
virtual photon as scalar particles. Here P = (M, 0) is the 4-momentum of a target (in the 
laboratory frame), P — p is the 4-momentum of the struck parton and p = {po,p) is the 
total 4-momentum of all other constituents (quarks and gluons) to which we refer as the 
“spectator”. The vertex F is a sum of all possible diagrams describing the virtual nucleon 
disintegration into quarks and gluons. Fig. ||. 

The scattering amplitude (the square block in Fig. |I]) satishes a Bethe-Salpeter equation, 
shown schematically in Fig. ^ 

T = V + VGqT, (2.1) 

where the Green’s function Go can be written as 

{p\GoiV)\p) = igs{p)go{V - p)5{p - p). (2.2) 

Here gs{p) is the spectator Green’s function, and go(V —p) is a Green’s function for a struck 
quark : 


4 


-p) 


1 


(2.3) 


{V — pY — m? + k 

where V is the total 4-momentum of the system {V = P + q for the kinematics in Fig. 
and m is the struck parton mass. The driving (interaction) term V in Eq. O is a sum 
of all irreducible diagrams which do not include the struck quark and the spectator in an 
intermediate state. Since quarks are conhned at large distances, this term is singular for 
{p ~ p'Y 0 instance, {p\V\p') ~ (p — p')~^ in the case of linear-rising conhnement). 

The corresponding Bethe-Salpeter equation for the vertex function T is obtained from 
Eq. (|2.1|) by taking the limit ^ = M^. Since the amplitude T factorizes near the 

nucleon pole 


T(P,p,p') 


r(p — p,p)r(p — p',p') 

P2 - M2 


(2.4) 


one hnds the following equation 


r = VGa(P)T, 


(2.5) 


which is shown schematically in Fig. 

By introducing the interacting (full) Green’s function 

G = Gq + GqV Go + ■ ■ ■ = —— (2.6) 

(jQ — V 

we can represent the structure function W{Q^, u) as 

W(Q\i^) = ilm /4(P,p)(p|G(P + ,)|p')4(P,p')P)?V, (2.7) 

71 J [ZTip 

where 


p) = 7 0 ^^ = 9oiP - p)r(P - p,p) 


(P — p)2 — 

is the relativistic bound state wave function. 


The lA of the structure function W corresponds to G —> Go in Eq. (|2.7|) . 
imation can be applied if the contribution from the interaction V is small. 


(2.8) 


This approx- 
However, this 
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may not be the case, in particular because of the infra-red singularity in V. The singularity 
can generate large corrections from the higher order terms of the expansion (|2.6|) . Yet, it 
does not imply that the correction terms would remain large in any other expansions of the 
Green’s function G. Indeed, the bound state wave function is generated by the same driving 
term V, Eq. (|2.5| ), so that the singular part of V can be compensated by some part of the 
mass and kinetic energy terms of in the full Green’s function G = [Gg ^ — V]~^, when the 
latter is substituted into Eq. (^T|) . In fact, one hnds from Eqs. (|2.5|) , the following 
relation 


|r-G„-‘(P)]g,|<J.)=0, 


(2,9) 


which in nonrelativistic limit corresponds to cancellations between the binding potential and 
the kinetic and binding energy terms in the Schroedinger equation. 

Eq. ( p.9|) suggests to expand the total Green’s function in (^?7|) in powers of the operator 


h = h{P) = V — Gq (P) instead of the lA expansion (|2.6|) , in powers of V. We thus obtain 


G = G + GhG + GhGhG + 


where G~^ = G~^{V,P) = G~^{V) - G~^{P). It can be rewritten explicitly as 


( 2 . 10 ) 


(p|G-‘(p,p)|p') = {p\v(P) - r(p)ip') +s.-‘(p)[(P -p)" - {P-pf]Hp-p') (2.11) 

One notes that in a general, the driving term V is not local, and therefore it may depend 
on the total 4-momentum V of the whole system (Fig. However, for a local interaction, 
{p\V{V)\p') =V{p — p'), so that the terms containing V in Eq. (1^.11) cancel. Then 


{p\G{'P,P)\p'} = i g,(p)!l(P - p, P - p)S{p - p'), 


( 2 . 12 ) 


where 


g{P - p, P - p) = 


(2.13) 


(P — py — (P — p)2 -f ie 
is a modihed, quasi-free. Green’s function of the struck parton with 4-momentum V — p. 


Unlike the Green’s function ^fg, Eq. (p(^) , the modihed Green’s function g, as well as G, 


6 











depends on the target 4-nionientuni P, which is related to the entire interacting system. 
The appearance of such an additional parameter is not surprising, since g has been designed 
to approximate the interacting Green’s function. Note that the pole in g does not appear at 
the quark mass (P — as in the free Green’s function go, but at an off-shell point, 

(p-p)2 = (P-p)2. 

Substituting Eqs. (|2.10|) , (p.l2|) into Eq. (|^ ) we obtain the structure function W in an 
expansion in powers of G = gsg for V = P + q. Gonsider the limit ^ oo and x = Q^/2mz/ 
= const. Then u = Q‘^/2Mx ~ Q^, I'/lql —>■ 1, and 


g + q - p, P - p) = {P - p + qY - {P - pf P ie 

= 2(M - po)u + 2p ■ q - + ie ^ 


(2.14) 


Therefore Eq. ( 2.10 ) represents in fact an expansion in powers of 1/Q^, and 


F = uW = -Im($|G + GhG + ■ ■ ■ |<1)) = Pq + 7 ^ + + 

TT 


( 2 . 16 ) 


Each term of this expansion can be represented by a modihed Feynman diagram. Take for 
instance the hrst two terms 


Pn = —Im / i 


. d^p 




TT 


(27r)"‘ {P — p + qY — {P — pY -|- ie 


(2.16a) 


Pi f .PpPp' 

= —Im / t- 


^iP,p)9s{p)hiP,p,p')gs{p')^iP,p') 


(2.16b) 


Q'^ 7T J (27r)® [(P — p-I-g)2 — (P — p)2 + ie][(P — p'-f g)2 — (P — p')^-I-ie] 

One can easily see that these terms correspond to the two diagrams in Fig. ^ where the 
Feynman propagator (fc^ — m^)“^ for the struck parton with 4-momentum k is replaced by 
-{k- qYY^. We mark it with 

Let us consider the hrst order term. Pi, Eq. ( |2.16b|) , which involves the interaction V. 
One gets from Eqs.(^T^, ( |2.12|) 

h99sj^) = [h,g\gs\^) = [V,g]gs\^) 

{^\g,gh= {^\gs[g,h] = -{^\gs[V,g] 
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(2.17a) 

(2.17b) 











Therefore Eq. (|2.16b|) can be rewritten as 


= —0“ Im [[g,r],g] 9 ,|<*) 


( 2 . 18 ) 


where the interaction enters throngh the double commutator of V and the Green’s function 
g. Using Eqs. (|2.13|) , (|2.14|) we can write Eq. (|2.18|) explicitly as 


jFi = Im 


.d^pdy 


^{P.p)9s{p) -p')fV{p- p')gs{p')^{P,p') 


{2tiY [2(M - pq)v + 2p - q-Q'^ + ief[2{M - p'q)u + 2p' ■ q-Q'^ + ie\‘^ 

(2.19) 


Here we hnd the factor {q{p — p'))^ in front of V, which removes the infra-red singularity, 
V ~ l/(p — p'y, in In fact, this factor reduces the contribution of to the structure 
function even for a non-singular interaction, provided that the scattering amplitude of a high- 
momentum quark peaks in the forward direction. One can show that the same procedure 
removes the infra-red singularity also in higher orders in the expansion (|2.15 ) for F. A 
direct evaluation of Fi for the case of a heavy spectator is given in Sec. 4, where we 
explicitly demonstrate that this term is small compared to Fq, even at moderate Q^. Now 
we concentrate on the hrst term, Fq, in the expansion (p.l5 ). 


III. THE LEADING TERM 


A. New scaling variable. 


Since our expansion (|2.15|) minimizes the hrst order correction term, the structure func¬ 


tion F can be well approximated by the hrst term, Fq, Eq. (|2.16aj) , which corresponds 
to the hrst graph in Fig. Actually, this graph represents an inhnite sum of diagrams 
corresponding to gluon and quark emission. Fig. |^, which results in logarithmic corrections 
to Bjorken scaling and has been studied in great detail [Q. Our treatment of these processes 
is not diherent from the standard approaches except for the modihed propagator g, marked 
by in Figs. |^, |^. Let us analyze the consequences of this modihcation. As an example 
we consider the two diagrams shown in Fig. The hrst graph corresponds to a process 
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with no gluon emission. For simplicity we take the spectator quarks to be a diquark of mass 
rris, thus 


9s{p) = 


— ml + ie 

Substituting Eq. (|0|) into Eq. ( |2.16a|) one hnds 


(3,1) 


= V 


- p + qf - {P - pjf], 


(3.2) 


where is a component of the proton wave function, Eq. (|2.8|) , with the vertex F replaced 
by F^°\ corresponding to nucleon disintegration with no gluon emission (Fig. ^). Integrating 
over po and neglecting the contribution from negative energy states (pair production) we 
get from the hrst ^-function: po = Ep = {ml + Then using (|2.14|) we can rewrite Eq. 

O as 




d^p 

(^ 


\{M - Ep)v + pq- QV2 


pdp 


i'>'"’(ipi)P 




(3.3) 


where = *h(T’, Ep^p)/2Ey‘^. The variable y is the minimal momentum of the struck 

quark, —Pmin = VQ/IqI^ obtained from the equation 


2 (m - \Jp + mA V - 2p|<j| - ^ 0. 


(3.4) 


The latter corresponds to energy conservation, given by the h-function in Eq. (|3.3|) , when 
the struck quark is equally off-mass shell before and after the virtual photon absorption. 
Solving Eq. (|^) and using i> = Q^/2Mx, we obtain 

(1 — x)‘^ — {ms/MY 


y{x,Q )/M = 


(3.5) 


Y^(l — xY + ^m\x^ -1- y^(l — xY + 4:M‘^x‘^{l — xY 

It follows from Eq. ( p.3| ) that {\q\/i^)T^'’ = f{y) is a function of the scaling variable y only. 
In the nonrelativistic limit and for zero binding energy, this variable coincides with the West 
scaling variable y, Eq. (El). 
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Eq. (^) for JF® can be rewritten in terms of light-cone variables: p = (p+,p_,p_|_), 
where p± = po ^ pz- Here the negative ^-axis has been chosen along the virtual photon 
direction, so that q± = =F IqI, and = 0. Introducing the light-cone fractions z = 
(F+ — p+)/P+ and ^ = —g+/P+, where P± = M is the target light-cone momentum, and 
integrating over p_, one obtains 


= Z/ 


(fp±dz 

2(27r)3 


6 {M‘^ + Mq_ 


m? + p]_ ml + p \' 
1-^ . 


(3.6) 


where Mg_ = Q'^H and 


= (P — pY = zM^ — 


2 zml + p\ 


(3.7) 


The expression ( ^.6| ) for clearly corresponds to the hrst diagram in Fig. calculated 
according the rules of light-cone perturbation theory, where the struck quark mass equals 
fh. Similar to the previous calculations one hnds after integration over and z in (p.6|) 
that the structure function depends only on a single (scaling) variable x. The latter is the 
value of the light-cone fraction z, corresponding to the zero of the ^-function argument in 
Eq. (|3.6|) for p_L = 0. One hnds 

X + \Jl + 4:M‘^x’^/Q‘^ — — xY + 


X = 


(3.8) 


1 + + 4M2a;7Q2 

Since x corresponds to the light-cone fraction of the off-shell struck quark with 4-momentum 
P — p (Fig. where p = {\Jml + —^9/1^1), one can easily verify that x is related to y 
by 


X = 1 — 


mf + + y 


(3,9) 


It is interesting to compare x with the Nachtmann |p scaling variable, ^ 

2x 


e = 


(3.10) 


1 + p+AhPx^lQ’- 

which corresponds to the light-cone fraction of the on-mass-sheh struck parton of zero mass 
I^ . The denominator of this expression takes into account the kinematical target mass 
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corresctions to the x-scaling. Similar target mass effects are accounted for in Eq. 
for the scaling variable x. It is not surprising, since no assumptions have been made for 
the value of the target mass M in Eqs. ( |3.3| ), ( |3.6| ). However, x includes also dynamical 
corrections to the x-scaling, which are not present in the variable These are taken into 
account by the last term in the numerator of Eq. ( |3.3|) . It can be seen in a more pronounced 
way if we consider large limit. Then using Eqs. (|^) , (|3.1(]|) one easily finds that x and 
^-variables are related by 






2 2 
m^x 


(3.11) 


Thus X —^ X for —>• cx). However, the term m^x^/(1 — x)Q‘^ makes x and ^ be quite 

different for finite in particular for high x. As an example we plot the variables ^ and 
X as functions of for x = 0.75 in Fig. H, for a spectator mass = M. One finds that x 
approaches the Bjorken variable x (the dotted line) much more slowly than the Nachtmann 
variable (^, and the difference between x and ^ is appreciable even for rather large Q^. For 
small X, however, the variable x is close to ^ or x, unless the spectator mass is not very 
large. 


B. Gluon radiation. 


Consider the second diagram in Fig. which describes gluon emission. Its contribution 
to the structure function can be written as 


•7r(l) 

0 



(27r)8 




- - p'f5[{P -p + qf 


{P-Pf] 

(3.12) 


where is the component of the proton wave function corresponding to the vertex T*^^) 
in Fig. |[ It is related to by 


|$d)(p. 


p,p-p)\ = 


47ra,|«l>W(P,p0p 

{p — p'Y — im? 


(3.13) 
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where a* is the QCD running coupling constant. After integration over the po and p'g we 
can rewrite in the form of Eq. (|3.3|): 


1^1 -r(l) _ 




(3.14) 


Here y{ms) is given by Eq. (3.5), where irtg is now the invariant mass of the diquark-gluon 
system: = p^. In terms of light-cone variables (see Fig. §) it can be written as 


2 2 , n ^{P'±-PA.f , / ! 

"i. = + (1 - z) —3—:— + - z)- 


Z' — Z 


l-z' 


/2 2 


(3.15) 


where rud is the diquark mass. Notice that 1 > z' > z. Using Eq. ( p.l5|) we rewrite Eq. 
( p.l2| ) in terms of light-cone variables 


d"’ = 


({‘pxdzd^p'^dz' ^ _m‘ + pi ml+pY 


4(27r)® {1 — z'){z'— z){z — ^) 




l-z 


(3.16) 


where = {P — p)^, Eq. 

Eq. (|3.16|) describes the lowest order gluon emission contribution to the structure func¬ 
tion, JFq. The total contribution of gluon emission to Pq corresponds to the sum of all ladder 
diagrams. It can be done using the same procedure as for instance in |lTI| , p!2| . Our modihca- 
tion consists only of replacement of the struck quark mass m by rh, Eq. (0. Eventually, 
it will modify the evolution equation []^ by the replacement of the Bjorken scaling variable 
X by X, Eq. (p.8|) , for rus given by Eq. (|3.15|) . 


C. Approximation. 

A considerable simplihcation can be achieved if we approximate as an effective spec¬ 
tator mass depends only on external momenta. Since z —^ x ~ x and z' ~ z, one gets from 
Eq. (|3T5|) 

~ 0(x, Q^)(l — x) (3.17) 


12 






















In this case the off-shell mass fh of the struck quark, Eq. (|3.7|) , becomes independent of 
z',p'j_. This allows us to integrate over z',p\ in Eq. ( |3.16| ), thus reducing it to the form of 
Eq. (|3.6|) (or ( p.3|) ). Finally one can sum all the ladder diagrams, such that (p±, z)\^ 
in Eq. ( p.6|) is replaced by u{p^,z,Q‘^), which is a probability to hnd a struck quark with 
momentum p inside the nucleon, and a spectator with any number of gluons. The latter 
gives rise to logarithmic corrections to scaling. 

Eq. ( p.l7|) for the invariant spectator mass looks quite appealing apart from its relation 
to Eq. ( 3.151) . Indeed, X = 1 corresponds to elastic scattering, when no gluons are emitted. 
Therefore in this case the spectator is represented by a diquark. When x decreases, gluons 
are emitted and ml increases oc (1—x). The coefficient C(x, in Eq. ( p.l7| ) determines the 
rate of increase of the spectator mass with and x. It can be found self-consistently from 
the evolution equation. However, when x ~ 1, one can take C'(x,(5^) — (7(1, ~ const, 

because of Q^-dependence of the spectator mass is less important than its x-dependence 
near the elastic threshold. Let us roughly estimate the value of C by using the Weizsacker- 
Williams or “equivalent photon” approximation, utilized in Ref. [Q for derivation of the 
evolution equation. One hnds from []^ that the light-cone fraction of the “equivalent” gluon, 
z — z\ (Fig. 1^) is of order asln((5^/Qo) region of large x. However, is also 


about the same order of magnitude. Then, as follows from Eq. (3.15|) , C r'U (b'x 
so that one could expect to hnd C on the scale of (GeV)^. 


IV. CORRECTION TERM 


Consider the hrst correction Eq. (|2.15|) , to the leading term Tq. In order to 

simplify the evaluation of Eq. (p.l9|) , we take the large spectator mass limit: S> |p|, 

so that Ps S{Po — In the same limit one gets from Eq. ( |3.4| ) 


y^y = 


+ 


m*v 


m\ 1^1 

where m* = M — m^, and also the Green’s function g, Eq. (3.14|), reads 


(4.1) 
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9 = 


‘2\q\{Pz + y + ie)' 


(4.2) 


Then the first order correction term Eq. (|2.19|) , becomes 

1 


JFi = —Ai Im 
where q = q/\q\, and 


27r J (27r)8 (p, + y + + i/ + !e)F(P ^ 


(4.3) 


Ai — 


Mx 


Mx 


q"^ Arusiy 2ms(l + 2m. 


, for Q‘^ ^ oo 


(4.4) 


It is convenient to evaluate the terms in (|2.15|) using coordinate representation. Sub¬ 
stituting 


4>{p) = / e*^'^0(r)(i^r, y{p — p')= I e 


/\ _ f J(P-p'yr 


v{r)d?r 


(4.5) 


into Eqs. (|3.3| ),( |4.3|) one obtains after some algebra the following expressions for the 
two first terms of the expansion (p.l5|) 

Fin 


J\q f 

jFg = — / dsexp{—iys) / — sq)(l){r) 

2ti J-oo J 


(4.6a) 


T\ = i 


.Fli roo 


/ oo r rs 

dsexp{—iys) d^r(j){r — sq)(j){r) / da[v{r — aq) — v{r)] (4.6b) 

-oo J Jo 

where Flo = ^/\q\ 1 for ^ oo. These expressions are of the same form (up to the 

coefficients Flj) as the first two terms of the nonrelativistic l/|g| expansion of the structure 
function H. One checks that the same corerspondence also holds for higher order terms 


jFj in the expansion (|2.15|) . Explicit analytical evaluations of the first three terms in the 
nonrelativistic l/|q| expansion §] for the harmonic oscillator and for the square well potential 
can be found in 0. 

Let us consider the case of the linear-rising potential, v{r) = yr. Notice that yji— aq\ — 
7 |r| ~ —^zajr for large r, and therefore the integrand in Eq. ([4.6b| ) does not increase with 
r. (It reflects the elimination of the infra-red singularity in the expansion ( |2.15| )). Using 
this result we perform the cr-integration in Eq. ([4.6b| ) and then relate T\ to the zero order 
term Tq by 
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^ ^2 /■- 


[ dsexp{-iys) [ d^r-(j){r - sq)(j){r) ~ 

J—oo J T ^Tf\ Oq 


(4.7) 


47r dy"^ J — OQ J r 2ro dy^' 

where rg is some average size of the system which determines the slope of the parton mo¬ 
mentum distribution. Approximating the structure function at small ?/ by a Gaussian, 

~ exp(—rQ?/^), and using Eq. (|4.7|) one obtains the estimate of the hrst correction 

term 


Q^:Foiy)^ 'Q2 


(4.8) 


where A = yro is of the order of the hrst excitation energy. It means that the correction 
term is indeed small, so that the approximation of the structure function by jFg should be 
valid also for non-asymptotic values of ^ few (GeV/c)^. 


V. COMPARISON WITH DATA 


In Fig. ^ we displayed the proton structure function from BGDMG [Q] and 

SLAG-MIT [0 experiments, as functions of for hxed x. The solid curves correspond to 
the 15-parameter £t to these data taken from |^: 


Fi{x,Q^) = A{x) 


lti(Qii/A2)J qG ’ 


(5.1) 


where Ql = 20 GeV^, A = 250 MeV, and 

n _ ^]‘2.562 4 0 0 1 5 ^ 

Aix) =-p—TT-j— 'Y' 0^(1 — xY] Bix) = 0.364 — 2.764x -I -; G(x) = 

x +0.0186’ ^ ’ tA 


where {a4=(0.4121, -0.518, 5.967, -10.197, 4.685), and {q}=(-1.179, 8.24, -36.36, 47.76). 

Let us treat this £t as the actual data and display it in Fig. 0 as a function of the scaling 
variable X, Eq. (|T8|): (the dashed lines). The spectator is taken 

to be a diquark of mass rUg = md = 850 MeV. This value of rud is taken from the bag model 
or the non-relativistic quark model, which estimate the scalar and vector diquark masses 
somewhere within the range of 600 to 1100 MeV [^,^. For a comparison we also displayed 
in Fig. 0, the structure function Eq. ( b.l|) , as a function of for hxed values 
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of X (the solid lines). The values of for x = const are taken within the limits of the data, 
i.e. x{x,Q^) < 0.75. Therefore the minimal values of for x=0.65 and 7=0.7 are smaller 
than the corresponding values of for x=0.65 and x=0.7. Since the values of x and 
for 7 = 0.75 are outside the kinematical region of the BCDMC and SLAC-MIT experiments 
{x{x, > 0.75), the case of x,x = 0.75 is not shown in Fig. 0 

One observes that Ff at constant 7 exhibits considerably weaker Q^-dependence than the 
same structure function taken at constant x. It implies that the essential part of the Bjorken 
scaling violations, usually attributed to the logarithmic terms, is in fact l/Q^-corrections, 
incorporated in the scaling variable 7. In order to assess what part of the x-scaling violation 
is accounted for by use of the scaling variable 7, we consider the following procedure. Let us 
take the proton structure function, F|’(7, Q^), at = 250 (GeV/c)^, which corresponds to 
the largest value of for the data sets in Fig. |^. In that region of the scaling variable 
7 is very close to x, and therefore F^i^x, 250) = Ff (x, 250) (see Fig. 1^). The corresponding 
“asymptotic” proton structure function /^(x) = F^x, 250), obtained from the £t (|5T|), is 
shown in Fig. O, together with three data points for = 250 (GeV/c)^. The dotted part 
corresponds to the same £t, Eq. (|5.1|), extended beyond the limit of the data (x > 0.75). 
The scaling of the structure function in the 7 variable, F|'(F, Q^) = /^(F), generates the 
Q^-dependence of the same structure structure function taken at constant x. 


Fnx,Q^) = rMx,Q% 


(5.2) 


Then the deviations of Eq. (^.2|) from the data would explicitly show of what part of the 
x-scaling violations is not incorporated in the variable 7. 

The results of this analysis are presented in Fig. |T^, with F^x, Q^), Eq. (^?2|) , given 
by dashed lines. As in the previous analysis (Fig. 0 the spectator has been taken to be a 
diquark, rUg = = 850 MeV and G = 0 in Eq. ( |3.17|) . One sees that the Q^-dependence 

of at X = 0.75 is well reproduced. The Q^-dependence of the other data sets is reproduced 
only partially, and the deviations from data increase for smaller x. However, the increase 
of spectator mass, rus, as 1 — x, Eq. (|3.17|) , can well influence the Q^-dependence even in 
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the region of small x. The evolution of the spectator mass with x and is given by the 
coefficient C in Eq. (|3.17|) . For large x this coefficient can be taken as a constant. It then 
may be extracted from data by requiring that Eq. reproduce two large x data sets, 

for instance x = 0.75 and x = 0.65. Since each data set is htted by adjusting the spectator 
mass Eq. (^.17| ) Exes also the parameter which is the value of diquark mass in the 
elastic limit (x = 1). Since in this limit the nucleon is not excited, one expects rud to be on 
the order of two constituent quark masses, i.e. 500 - 1000 MeV jlBHIHl. It appears however, 
that the two data sets cannot be fitted with such values of m^, but only with 0 and 

C ~ 3 (GeV)^. This corresponds to a spectator built out of very light quarks. 

Taking rud = 0 in Eq. (|3.17|) , we find for the scaling variable x, Eq. 

X + \Jl + — y^(l — xY + 4G(1 — xjx^/Q^ 


X = 


(5.3) 


1 + +4M2xVg2 

One gets from Eq. (|5.3|) that x = 1 for x = 1, so that the two scaling variables vary within 
the same limits. 

Proton structure function given by Eqs. ( |5.2|) , (|5.3|) for G = 3 is shown in Fig. 0by 
dashed lines. Rather good agreement with the data is observed even for x < 0.65, although 
one expects large logarithmic corrections in this region and also the variation of C should 
be taken into account when x is far from 1. It is therefore of greater interest to make a 
comparison for x > 0.75. We show in Fig. 0a the data for the proton structure function 
taken from new SLAG measurements in the threshold region for 7 < < 30 (GeV/c)^ |l9 


together with three high-statistics spectra for g^=5.9, 7.9, and 9.8 (GeV/c)^ from a previous 
SLAG experiment These data do not scale either in the variable x (Fig. p!4|a), or in the 
Nachtmann variable ^ |^. In contrast, excellent scaling (Fig. |T^) is observed when the 
data are plotted as a function of x, Eq. ( b.3|) , with G=3. It is even more remarkable that 
the structure function obtained from these data completely coincides with the asymptotic 
structure function /^(x) = T 2 ’(x, 250) taken from BGDMG experiment (Fig. |^) and shown 
by the solid line. This confirms the dominance of the quasi-free term JFq; Eq- ( |2.15|) , which 
is in line with our estimates in the previous Section. Indeed, the lowest correction term, jFi, 
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is proportional to y, Eq. (^4.8|) . Since y = Q ioi x = 1 and iris = 0, Eq. (|3.5|) , it follows that 
> 0 near the elastic threshold. 

The analysis of the large x data ||T^ , |20[] (Fig. |T^) allows us to extract the asymptotic 
structure function f^{x) = up to x ^ 0.95. We hnd that f^{x) is clearly below 

the dashed curve for x > 0.8, which is the £t (b.l|), extended outside the data, Fig. 0- 


It means that the fit (|5.1|) is not applicable in that region. It was shown by Drell, Yan, 
and West [0,^ that the threshold (x —> 1) behaviour of the asymptotic structure function 
is correlated with the large behaviour of the elastic form factor. Using quark counting 
arguments one obtains /^(x) ~ (1 — x)^. However additional QCD effects may effectively 
increase the value of the exponent [^. This predictions can be checked by a comparison 
with the structure function that we extracted from the data. We hnd that it well described 
by /^(x) = 1.5(1 — x)^’’^ for x ^ 0.75 (the dot-dashed curve in Fig. |T^). 

With the asymptotic structure function /^(x), extracted from the experiment |]I|,p!^,^ , 
we can predict the structure function at large values of x from moderate up to very high 
values of Q^. The results are shown in Fig. H- The dashed lines correspond to Eq. o. 
with the asymptotic structure function /^(x) given by the £t (^.1| ) for x < .78, = 250 

(GeV/c)^, and /^(x) = 1.5(1 — x)^’’^ for x > 0.78. We plot also a few available data points. 


taken from old SLAC-MIT measurements for x = 0.78, 0.82 and 0.86. Again our 
predictions are in full agreement with the data. Still a check of our predictions for higher 
would be of greater interest, since we predict a signihcant non-logarithmic Q^-dependence 
for large x. Fig. |^. In fact, as follows from our analysis, none of the large x data (x > 0.5) 
exhibit any substantial deviations from the x-scaling. Therefore, measurements for high 
would be extremely important in order to establish the importance of logarithmic corrections 
in the large x region. 
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VI. SUMMARY 


In this paper we concentrated on effects of conffnement in deep inelastic scattering. Using 
the framework of the Bethe-Salpeter equation, we found a new expansion of the structure 
function in powers of 1 /Q^ that is free of infra-red singularities and diminishes corrections 
to the zero-order term. The zero-order term describes scattering off a free off-mass-shell 
parton, which keeps the same off-shell mass in the final state. We evaluated corrections 
from higher order terms for the case of a linear-rising confining potential, and found them 
small even for rather low values of Q^. It allows us to consider the zero order (quasi-free) 
term as a good approximation for the structure function, valid in the entire region. 

By analyzing the quasi-free term we found it depends on the scaling variable x. This 
variable coincides with the Bjorken variable x in the limit oo. However, at finite 

the variables x and x are quite different: x — x ^ ~ x)Q‘^, where rris is the invariant 

mass of the spectator particles. It implies that l/Q"^ corrections to the x scaling would 
be much larger than those obtained in perturbative calculations, especially in the region of 
large x. These corrections could be very appreciable also at small x, since the spectator 
mass increases with [l—x) due to gluons emission. However the evaluation of the spectator 
mass at small x depends on knowledge of its Q^-dependence, which can be obtained by using 
the evolution equation. In this paper we limited our analysis to the large x region, where 
the Q^-dependence of the spectator mass is less important. 

Using simple arguments we showed that ml = + U(1 — x) for a; ^ 1. Here is the 

diquark mass, since there is no gluon emission at the elastic threshold. First we analysed the 
data for the proton structure function by taking for values between 500-1000 MeV, i.e., 
on the order of two constituent quark masses, and taking U = 0. Even though we found that 
scaling in x is certainly better than that in x, the scaling deviations are still considerable, 
especially when x approaches the elastic threshold. Next we included the gluon emission 
contribution to the spectator mass by taking U 7 ^ 0. However, instead of taking for the 
diquark mass, md, the values from constituent quark models, we considered the nucleon 


19 



structure function data as a source of information for the value of We found that all 
large x data {x > 0.5) display perfect scaling in the x-variable for C ~ 3 and ~ 0, which 
corresponds to very light quarks. 

Since x ^ x for cx), the perfect scaling in the x-variable allows us to arrive at the 

Bjorken limit already at moderate values of . Thus our analysis of the proton structure 
function near elastic threshold shows that F 2 ~ (1—for a; —1. This is different from the 
theoretical results based on simple quark counting arguments, which predict T 2 ~ (1 — a;)^. 

Untill now the measurements of the structure function for large values of a; (a; > 0.75) 
have not been extended beyond ~ 25 (GeV/c)^. The present large x data are in full 
agreement with our predictions based on a;-scahng, and thus do not display any noticeable 
effects of logarithmic terms. It would be very interesting to extend the measurements to 
higher values of Q^. The a;-scaling predicts a distinctive Q^-dependence in the structure 
function F{x, Q^) at large hxed values of the Bjorken variable x. Hence, the deviations from 
our predictions would establish the magnitude of the logarithmic scaling violations at large 

X. 
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FIGURES 


FIG. 1. Nucleon structure function given by the imaginary part of the forward Compton am¬ 
plitude. The first diagram is the Impulse Approximation, and the second one describes the Final 
State Interaction. The shaded area includes spectator particles (quarks and gluons). 

FIG. 2. Nucleon vertex function, which describes quark and gluon emission. Quarks are shown 
by solid lines and gluons by wavy lines. 

FIG. 3. Proton structure function from BCDMS [Q| and SLAC-MIT Q experiments. The solid 
curves correspond to a 15 parameter fit to these data, taken from [^. 

FIG. 4. Bethe-Salpeter equation for the operator T describing interactions between the struck 
parton and the spectator partons in the final state. 


FIG. 5. Bethe-Salpeter equation for the vertex F describing the relativistic bound state. 


FIG. 6. Diagrammatic representation of the first two terms of the expansion ( 2.15| ). Modified 
propagator of the struck parton is marked by . 


FIG. 7. Diagrammatic representation of the leading term. Quarks and gluons are shown by 
solid and wavy lines respectively. The modified propagators are marked by . 

FIG. 8. The same as in Fig. but using light-cone variables. The negative 2 ;-axis has been 
chosen along q. 


FIG. 9. Gomparison between the scaling variable x = x{x,Q‘^) and the Nachtmann variable 
C = ax,Q^), Eqs. O, for X = 0.75. The spectator mass equals the nucleon mass. The 

dotted line corresponds to x=0.75. 


FIG. 10. The proton structure function plotted as a function of for fixed values of the 
Bjorken variable x (solid lines) and the scaling variable x, Eq. (|3.8|) (dashed lines). 
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FIG. 11. Proton structure function in the asymptotic region, f^{x) = (x,250), given by the 

fit, Eq. O- Three data points are taken from Q . The dashed part of the curve lies in the region 
outside the data Q. 


FIG. 12. Q^-dependence of the structure function F^ix^Q'^), which corresponds to scaling in 


the variable x, Eq. (^). The spectator is taken to be a diquark of mass mg = md=850 MeV. 


FIG. 13. The same as in Fig. 12, but where the spectator mass is a function of x, Eq. (^.17 ) 
for md = 0 and C = 3 (GeV)^. 


FIG. 14. The structure function in the region of large x for 7 < < 30 (GeV/c)^ 


20 |, 


plotted (a) as a function of the scaling variable x and (b) as a function of the scaling variable 
x, Eq. (|5.3[) . Three high-statistics data sets |^] for Q^=5.9, 7.9, and 9.8 (GeV/c)^ are marked 


by “x”, and respectively. The solid curve and the dashed curves correspond to the 

asymptotic structure function, shown in Eig. |Tl|. The dot-dashed curve in (b) corresponds to 
F|(x,g2) = 1.5(1 -x)3-7. 


FIG. 15. Predictions for the structure function F 2 {x,Q‘^) in the region of large x. The data 


points @ correspond to x =0.78 (-I-), 0.82 (x), 0.86 (#). 
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